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Abstract: The prospects for making a direct measurement of the W production charge 
asymmetry at the LHC are discussed. A modification to the method used at the Tevatron is 
proposed for measurements at the LHC. The expected sensitivity for such a measurement 
to parton distribution functions is compared to that for a measurement of the lepton 
charge asymmetry. The direct measurement approach is found to be less useful for placing 
constraints on parton distribution functions at the LHC than a measurement of the lepton 
charge asymmetry. 
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1. Introduction 

The measurement of the W production charge asymmetry at hadron-hadron colliders pro- 
vides important information about parton distribution functions (PDFs). The simplest 
method is to measure the asymmetry in the pseudorapidity distribution of the charged 
leptons (the lepton asymmetry) arising from leptonic decay of the W^ bosons. Such mea- 
surements have been performed at the Tevatron, for both W — > ev and W —> [iv events, 
by both the CDF [l], |2| and DO ||, |4| collaborations and the data have been included in 
global fits of parton distributions |], |(| [/J . The asymmetry is mainly sensitive to valence 
quark distributions ||, providing complementary information to that obtained from mea- 
surements of structure functions in deep inelastic scattering [||, 10, 11, 12, [H|. Predictions 



for the asymmetry exist up to next-to-next-to- leading order in the strong coupling [[Tj] . 

The CDF collaboration has recently published a direct measurement of the W produc- 
tion charge asymmetry [15]. This measurement used a technique first proposed by Bodek 



et al. 1 16], in which the W boson rapidity is inferred event-by-event from the lepton four- 
momentum and the missing transverse momentum. This inference enables a measurement 
of the W asymmetry to be made at values of rapidity where the uncertainty on the pre- 
diction from theory is much larger than the uncertainty on the lepton asymmetry for any 
part of the pseudorapidity range. 

The inclusion of the CDF W asymmetry [15] in PDF fits has not been straightforward: 
the data set gives the second largest contribution per data point to the \ 2 °f the NNPDF 
fit Q. Attempts by the MSTW group [17] to include it in their fits have revealed tension 
with the electron asymmetry measured by DO ||. To investigate this the CDF collaboration 
have produced the electron asymmetry for the same dataset as used in the W asymmetry 
measurement with two interesting features revealed. Firstly this electron asymmetry agrees 
well with the DO electron asymmetry and secondly, while the W asymmetry agrees very well 
with the predictions made using the CTEQ6.6 [jig] PDF set, the CDF electron asymmetry 
disagrees with predictions made using the same PDF set. This apparent inconsistency may 
indicate biases which have not been accounted for in the method. 

In this paper the feasibility of measuring directly the W production charge asymmetry 
at the LHC is investigated for the first time using simulated data. The expected sensitivity 
to the PDFs for such a measurement is compared to that for a measurement of the lepton 
asymmetry. 

The paper is structured as follows: the simulated data samples used for this study are 
summarised in section ||; the procedure used to measure the W asymmetry is described 
in section El. Problems specific to the LHC environment are described and modifications 
to the reconstruction procedure proposed in sections || and ||. The performance of the 
modified method is studied in greater detail and the expected statistical and systematic 
uncertainties on the W asymmetry are compared to the expected statistical uncertainties 
on the lepton asymmetry for different luminosities in section [?]• Finally, the implications 
of this study for early measurements of W production charge asymmetries at the LHC are 
discussed in section ||. 

2. Simulation of W production 

The studies in this paper are based on samples generated using the Monte Carlo (MC) 



generators Pythia 8.120 |19| and Herwig++ 2.4.0 [g0| glj. The Pythia sample is a 
simulation of the process //' —> W, where / and /' are fermions, at leading order (LO) in 
the strong coupling. The Herwig+- 1- sample is a simulation of the W production process 
at next-to-leading order (NLO) in the strong coupling. As such, it represents the state- 
of-the-art in simulation of Drell-Yan vector boson production. The positive weight NLO 



matching scheme (POWHEG) |22| , 23] was used in the generation of the Herwig++ sample. 
This approach consistently combines the NLO calculation and parton shower simulation 
and was implemented in Herwig+- 1- by Hamilton et al. |24[| . 

The PDFs used for the generation and the size of the samples used are specified in the 
appropriate sections of the paper. No detector simulation is applied. The following fiducial 
cuts are used in the event selection in order simulate the detector acceptance: 



a cut on the transverse momentum of the lepton, Pj, > 25 GeV, 

a cut on the missing transverse momentum (the transverse momentum of the neu- 
trino) in the event fix > 25 GeV, 

a cut on the lepton pseudorapidity \r) l \ < 2.4. 



3. Analysis Technique 

The procedure used to extract the W production charge asymmetry as a function of the W 



boson rapidity, yw, follows closely that used by Bodek et al. [16]. It can be broken down 
into the following steps: 

1. Calculation of the W rapidity solutions: For each reconstructed event, the 
two possible W rapidity values (referred to as "solutions") are calculated from the 
missing transverse momentum in the event and the charged lepton momentum using 
a constraint on the W boson mass. 

2. Weighting of the W rapidity solutions: Event-by-event, the solutions are filled 
into separate y w + and y w - histograms, using weights calculated from MC input as 
described in section 01. If only one physical solution is found, the solution is filled 
with a weight of 1. Step 1 (calculation of the solutions) and step 2 (weighting of the 
solutions) are hereafter referred to as full kinematic W reconstruction. 

3. Acceptance corrections: The rapidities of the kinematically fully reconstructed 
W bosons are corrected in each bin of yw for the detector acceptance and also for 
biases of the full kinematic W reconstruction arising from the weighting procedure 
(as discussed in section |]). The acceptance corrections are calculated bin-by-bin as 
the ratio between the generated yy/± distributions after cuts and full kinematic W 
reconstruction and the MC generated W distributions before cuts: 



MC events after cuts and with full kinematic W reconstruction(bin n 

S(bini) 



^ m l ' MC events before cuts/ 



The acceptance correction is applied by multiplying 1/Acc( bini ) with the content of 
bin i of the yw distributions of the fully reconstructed W bosons. Step 3 yields the 
reconstructed yw distributions, extrapolated to the whole phase space. 

Compare reconstructed and MC-input yw distributions: If the reconstructed 
yw distributions and the MC input yw distributions are consistent with each other 
then the procedure stops (go to step 7). If they disagree, then the procedure must 
be iterated, with the MC input reweighted to describe the reconstructed distribution 
(go to step 5). 



5. Reweight input MC to reproduce the reconstructed yw distributions: The 

input MC yw distributions without any cuts applied and the yw distributions recon- 
structed in step 1-3 are compared. A reweighting factor is extracted bin-by-bin: 

_ experimentally determined y w (bin i} 
r[ lj " Input MCy^ e (bini) {d ' 2) 

Each event in the input MC is re-weighted with the event weight r (y w ue (MC input)). 

6. Recalculation of weights: The weights used in step 2 and the acceptance correc- 
tions (step 3) are recalculated from the reweighted MC input and the steps 1-4 are 
repeated until convergence. This iteration reduces the dependence on the MC input. 

7. Measurement of the W asymmetry: When the reconstructed and modified input 
MC yw distributions are in agreement within their statistical uncertainty, the W 
asymmetry is considered to be determined. 



4. Pull kinematic W reconstruction: weighting of the W rapidity solutions 

In the weighting procedure |lgj (step 2 from section y) the twofold ambiguity for p u z and 
thus for yw is resolved statistically with the help of MC predictions. The possible rapidity 
solutions si and S2 are weighted with their respective probabilities (Pi and P%) and the 
total probability is normalised to unity, V = P1 + P2 = 1- The probabilities of each solution 
occuring are derived, based on: 

1. Expected cross sections da/dyw'- The cross sections da/dyw have a distinct 
behaviour as a function of yw- The two different rapidity solutions can be compared 
as to which is more probable. The expected cross sections da/dyw predicted by 
Monte Carlo are used as probability density functions to determine this. 

2. Expected lepton decay angle and anti-quark/quark ratios: The lepton decay 
angle in the W rest frame, cos6*, follows a (1 ± cos9*) 2 distribution, which can be 
used as another probability density function to weight the two rapidity solutions. 
The sign of the (1 ± cos6*) 2 distribution depends on the helicity of the lepton and 
the helicity of the incoming parton with highest-x, hereafter referred to as "higher- x 
(anti-)quarks". Therefore in the cos#* weighting, the probability density function is 
built from the contributions of events with higher-x quarks (l=Fcos 0*) 2 and of events 
with higher-x anti-quarks (1 ± cos#*) 2 . The probability density function for the 
weighting is constructed as the sum of the two decay distributions. The contribution 
of the higher-x anti-quark events is adjusted using the parameter -, which is defined 



as 



q w number of higher- x anti-quarks w , . 

-{PT,yw) = r f-tt-t ; — (Pt ,Vw) (4.1) 

q number of higher- x quarks 



The ratio | is a function of pYp and yw, it is determined from Monte Carlo simulations 
and therefore also depends on PDFs. Using |, one can build the total probability 
density function from the sum of the two possible angular distributions for the LHC 



P(cos 6*) = (1 T cos 6*) 2 + -(p^,y W )0- ± cos 6*f (4.2) 



and similarly for the Tevatron [16]. 



The probability density function for the cross section and the decay angle can be 
combined to calculate weighting factors, w\ and u>2, for s\ and S2, depending on the charge 
of the W ± : 

P(cos 6*± w± ) x (da/dy w ± ) 

Eti (P(cos^ ±>w± ) x (da/dy w ±)) 

In a reconstructed event each of the solutions 1 and 2 is weighted by w\ and u>2 
respectively and filled with that weight into histograms of yw±. 

5. Performance of the kinematic W reconstruction 

The performance of the weighting and reconstruction of the W, steps 1 and 2 in the 
iterative procedure, was studied. This study was performed using 8xl0 6 Pythia MC 
events generated with the MSTW08 NLO PDF set §] for Tevatron (pp, y/s = 1.96 TeV) 
and LHC conditions (pp, y/s = 14 TeV). The analysis technique was tested using the same 
MC sample as MC input and as 'pseudo-data'. The transverse components of the neutrino's 
momentum, p v x and p v y were taken directly from the MC. No acceptance corrections were 
applied for this part of the study. 

5.1 Performance at the Tevatron 

Figure [j] shows the weight assigned to the rapidity solution calculated in step 1 that was 
closest to the true rapidity. The figure shows the normalised distribution only for W~ 
events but is very similar for W + events. The distribution exhibits a peak at 1 and lies 
mostly above 0.5, indicating that in the weighting procedure the solution closest to the 
true yw value is given the larger weight. 

Figure [2] a) shows the asymmetry distributions for the Tevatron. Shown are the true 
W asymmetry distributions from the MC with and without cuts. Also depicted is the W 
asymmetry that arises when only the W solution closest to the true yw value is used and the 
W asymmetry based on the full kinematic reconstruction. All asymmetries are in agreement 
up to yw = 1.5, beyond that there is some disagreement between the reconstructed and 
the true W asymmetry. This can be better seen in figure b), which shows the ratio 
of the fully reconstructed W asymmetry (without acceptance corrections) and the true 
asymmetry. The difference between the reconstructed and the true asymmetry in the 
region defined by the acceptance cuts is only around 20%. The true asymmetry before 
cuts is at most 80% to 100% larger than the reconstructed asymmetry. Hence it can be 
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Figure 1: Normalised distributions of weights given to the rapidity solution that is closest to true 
W~ rapidity in the Tevatron environment. 
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Figure 2: This figure shows the asymmetry as evaluated at various stages of the iterative weighting 
procedure (a). In b) the ratios of reconstructed asymmetry and true asymmetries (before and after 
applying cuts) are depicted. 



seen that acceptance corrections correct not only for genuine detector acceptance effects 
but also for biases of the full kinematic reconstruction. However, the former are larger and 
more important than the latter. 

5.2 Performance at the LHC 

It is not immediately obvious that it is possible to transfer the full kinematic reconstruc- 
tion technique to the LHC because of crucial differences between the Tevatron and LHC 
environments. The Tevatron is a pp collider operating at centre-of-mass energy y/s =1.96 
TeV, while the LHC is due to collide pp beams at y/s up to 14 TeV. These differences pose 
a major challenge when applying the weighting method to the LHC environment. Collid- 
ing pp beams, will break the mirror symmetry between W + and W~. The asymmetry in 
pp collisions is the same at positive and negative rapidities. At the LHC W + bosons are 
produced almost twice as often as W~ and generally with a larger boost along the Z axis. 
Hence the performance of the method will differ in the two environments. 



The full kinematic reconstruction and the weights used were studied using a Pythia 
MC with the MSTW08 NLO PDF set. As already described above, the two solutions 
calculated in step 1, were weighted in step 2 using two basic considerations. 

• Shape of cross section as function of yw' At the LHC the cross section for 
W production is much flatter and expands over a much larger range of yw than at 
the Tevatron. This is shown in figure [|. There is no peak in the yw distributions 
for central rapidities and the cross section begins to fall only above yw > 1-5 — 2.0. 
Therefore, if both reconstructed solutions are within —1.5 < yw < 1-5, both are 
equally probable. In particular, the reconstruction of W + rapidities suffers from this 
fact, since the longitudinal boost of W + bosons is generally larger than the boost of 
W~ , and therefore the rapidity distribution extends to higher yw- 

• Ratios of leading q versus leading q: A further problem is the contribution of Ws 
produced with a higher-rc anti-quark. If the higher- x parton participating in the Drell- 
Yan W production is an anti-quark, the W is no longer produced with a boost parallel 
to the incoming quark, but antiparallel to the incoming quark. This introduces a 
change of sign in the expected cos 9* decay angle. In the weighting procedure this 
is accounted for by constructing the cos 9* weight such that it combines a weight 
for leading quark and leading anti-quark W production according to their relative 
contributions by using the ratio q/q as a function of pjf and yw, as in equation L3. 
At Tevatron centre of mass energies, leading quark production is indeed the most 
probable process, with q/q=0.25 at most. At the LHC however, the ratio q/q even 
extends to values as high as 1 as shown in figure y. This introduces a strong ambiguity 
to the cos 9* weighting. This effect is particularly problematic for W~ production 
because the ratio of anti-quark induced processes is higher than for W + production. 
The effect is also shown in figure |], where the cos 9* distributions of W~ (left) and W + 
(right) bosons are shown for W production at the LHC. For the W~ distributions, 
the contribution of higher- x anti-quarks is especially large and two solutions with 
cos^*' 1 ~ — cos6>*' 2 cannot be distinguished. Solutions with cos6>*' 2 ~ — cos^*' 1 or 
cos9*' 2 ~ cos^*' 1 are in fact most common and in these cases both solutions will be 
weighted with factors of about 0.5. 

Figure |B| displays the weight given to the rapidity solution that was closest to the 
correct rapidity at the LHC. While at the Tevatron, this distribution peaks at 1 and is 
mostly larger than 0.5 (as shown in fig. ||), at the LHC most weights are around 0.5. For 
W~ production the weights given to the right solution are on average at least marginally 
larger than 0.5, however for W + rapidities it is narrow and symmetric around 0.5, thus 
reducing the impact of the weighting procedure. The larger spread of W~ weights comes 
from the different expected da/dyw rather than decay angle effects. 

Again it is instructive to evaluate the asymmetry distribution at the various stages of 
the iterative weighting procedure. This is done in figure y a). Additionally, figure y b) 
depicts the ratio of the fully reconstructed asymmetry distribution to the true asymmetry 
distributions. The figure shows that in the central region, the corrections for intrinsic 
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Figure 4: cosf5* distributions at the LHC. The total of the higher- a; quark (dotted line) and 
higher-x anti-quark (dashed line) distributions is shown. 



systematic biases of the full kinematic reconstruction dominates, while in the forward region 
genuine acceptance corrections dominate. Comparing this figure to the equivalent plots for 
the Tevatron (figure b) reveals that the acceptance corrections are more important at the 
LHC than at the Tevatron. While at the Tevatron the ratio of the reconstructed asymmetry 
and the true asymmetry with cuts is between 0.8-1.1, the same ratio is 0.5-3.0 at the LHC. 
Equally the ratio of the reconstructed asymmetry and the true asymmetry without cuts is, 
at 1.5-3.0, larger at the LHC than at the the Tevatron (0.8-1.8). 
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Figure 5: The weight given to reconstructed rapidity solution closest to the true value of yw for 
W~ and W+ at the LHC. 
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Figure 6: a)the asymmetry as evaluated at various stages of the iterative weighting procedure, b) 
the ratios of reconstructed asymmetry and true asymmetries (with and without cuts). 



6. Modification of the analysis technique for the LHC 



At the LHC the full W kinematic reconstruction suffers because in this kinematic region it 
is difficult to favour one rapidity solution over another. Even with the best possible weights, 
it would not possible to restore the true asymmetry perfectly. At both the Tevatron and the 
LHC, acceptance corrections compensate for intrinsic kinematic biases in the full kinematic 
reconstruction. As a result of these features a new simplified analysis technique for the 
LHC can be proposed: in this new simplified approach, the construction of weights for 
the full kinematic reconstruction is bypassed and instead each solution is filled with a 
weight of 0.5, equivalent to randomly choosing either solution. This approach relies then 
solely on the acceptance correction, which is then also the only step dependent on MC 
input. This simplifies the determination of the systematic errors of the method. The 
acceptance corrections should be applied in bins of yw and pjf. The iteration procedure 
with reweighting of the MC input is still necessary to remove the MC dependence. This 
approach is hereafter referred to as the "new scheme". 



7. Performance of the modified analysis technique at the LHC 

The expected systematic uncertainties of a measurement in the new scheme were evaluated 
in comparison to the uncertainty on the predicted VF-asymmetry and in comparison to 
the lepton asymmetry. Here, only the expected statistical errors and the inherent uncer- 
tainty of the modified W measurement method are considered, other potential sources of 
experimental uncertainty are not taken into account. 

7.1 Uncertainties on the acceptance corrections due to statistical fluctuations 

It is not trivial to propagate the uncertainties from statistical fluctations in the data sample 
and the acceptance corrections to a total statistical uncertainty on the W asymmetry. The 
method reshuffles the number of events in the bins and might amplify statistical fluctuations 
as was also observed for the original implementation of the method pq] . 

In order to estimate the statistical uncertainties on the W asymmetry, 450 toy ex- 
periments were carried out (375 in the case where the toy experiments contained 5 x 10 6 
events each). Two Pythia samples generated with the CTEQ66 PDF, y/s = 14TeV were 
used. The MC input sample consisted of 18xl0 6 events and was used to calculate the 
acceptance corrections. 450 sub samples of 1 x 10 5 , 5 x 10 5 , 1 x 10 6 and 5 x 10 6 events 
each were drawn randomly from the pseudo-data sample of 75xl0 6 events, corresponding 
to integrated luminosities of 5.4, 27, 54 and 270 pb _1 respectively. These sub-samples were 
not exclusive, no attempt was made to exclude events already used in one toy data sample 
from other toy samples of the same luminosity. 

For each of the toy "data" samples, the W rapidity solutions were calculated for each 
event and filled with a weight of 0.5. Next the acceptance corrections derived from the 
MC input sample were applied. The spread of the measured asymmetries in each bin of 
yw was then interpreted as the statistical uncertainty on the measured W asymmetry. 
Figure [?] a) shows this spread for the toy 'data' samples. The entries belonging to different 
bins in yw are shown in different grey tones. Figure [7] b) depicts the W asymmetry 
with the uncertainty extracted in the toy experiments as bands for the different integrated 
luminosities. Figure c) shows the relative uncertainties on the W asymmetry as a function 
of y w . 

7.2 Uncertainites on the acceptance corrections from PDFs 

Uncertainties on the acceptance correction arising from the chosen MC yyy± input distri- 
butions differing from those in data due to the underlying PDF distributions were also 
considered. In this study the CTEQ66 and MSTW08 PDF distributions were used. To 
estimate the uncertainty the PDF error sets based on the standard Hessian technique § 
were used. 

Figure M shows the measured asymmetries and the uncertainty on them arising from 
PDF uncertainties for the CTEQ66 and the MSTW08 PDF sets, both shown as a ratio to 
the central value of CTEQ66. Around the PDF predictions and their uncertainties, a dot- 
dashed envelope is drawn, which is taken to give the PDF uncertainty on the asymmetry 
measurement due to the acceptance corrections. This error band is also centred around 
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Figure 8: The uncertainties on the measured W asymmetry from statistics and from the PDF 
uncertainties in two different bins of p^ . 

the CTEQ66 value and drawn as a grey band. Figure ||] shows on the left the relative PDF 
uncertainty on the asymmetry for 1.15 < p^ < 2.47 GeV and for 11.3 < p^ < 24.1 GeV on 
the right. For these and also the other low and mid-p^p bins the uncertainties are between 
5-12%, for larger values of pjf, the uncertainty tends to decrease. 

7.3 Experimental uncertainty due to the detector 

In order to estimate the effect of the reconstruction of the W events in the detector, the 
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four-vectors of the generator level particles were smeared to approximate the effect of 
differing MC and data resolutions. This was done by applying a Gaussian resolution factor 
to the p x and p y variables of the neutrino and to the p x , p y and p z variables of the lepton. 
These resolution factors were drawn randomly and independently for each of the variables 
from a Gaussian distribution centered around one with a width of 1% (lepton) and 5% 
(neutrino) which represent realistic uncertainties on the resolutions at the LHC ^6|. This 
was done 500 times for the same ~ 16 x 10 6 Pythia events. The reconstructed asymmetry 
was corrected with the nominal acceptance corrections derived from the unsmeared sample. 
In a sizeable number of events, as a result of the smearing, no physical solution could be 
obtained for yw- In these cases, the neutrino pt is scaled down in steps of 0.1 GeV until a 
valid solution is found, following the prescription of the CDF measurement, where the $t 



was assumed to be overestimated and scaled down [15|. 

Due to the fact that the pt spectrum of leptons is correlated with 77, there is a slight 
bias in the mean of the smeared lepton asymmetry distributions with regard to the nominal 
(unsmeared) distribution. A similar but larger bias can also be observed in the measured 
W asymmetry distribution. When calculating the experimental uncertainty, the bias of the 
mean of the 500 smeared distributions with regard to the nominal distribution is added in 
quadrature to the RMS of the 500 smeared distributions for both the lepton and the W 
asymmetry. 

7.4 Application of the modified scheme to pseudo-data 

In the following the weighting procedure of the modified scheme is investigated for the 
LHC environment. Two cases were studied: 

1. Different PDFs: A Pythia sample generated using the MSTW08 PDF (8 x 10 6 
events) was used as pseudo-data and a Pythia sample generated using the CTEQ66 
PDF (75 x 10 6 events) was used as MC input. 5 iteration steps were used. 

2. Higher orders: An Herwig+- (- W production sample generated using the MSTW08 
PDF (1 x 10 6 events) was used as pseudo-data and a Pythia sample generated using 
the MSTW08 PDF (18 x 10 6 events) was used as MC input. 5 iteration steps were 
used. 

Figure || a) compares the true W asymmetries for the CTEQ66 MC input and the 
MSTW08 pseudo-data to the reconstructed W asymmetries as obtained in the iterations 
steps 1-5. The reconstructed W asymmetries from all of the iteration steps agree better 
with the MSTW08 pseudo-data than with the CTEQ66 MC input, showing that the method 
still works. This can also be seen in figure || b) which shows the ratio of the reconstructed 
W asymmetries at each iteration step to the true VF-asymmetry of the MSTW08 pseudo- 
data. The first iterations agree within 7%, for the other iterations larger deviations are 
visible. This is an effect of statistical fluctuations, the iteration procedure should be sensibly 
stopped after the first or second iteration for a measurement. This can be seen in figure ^ c), 
which shows the ratio of the reconstructed W asymmetries of the iteration steps to the W 
asymmetry of the specific previous iteration step or the input MC asymmetry in the case of 



12 



. 0.22 

0.2 

0.18 

' 0.16 

0.14 

0.12 

0.1 

0.08 
0.06 
0.04 
0.02 



MSTW08 Original asymmetry (pseudo-da' 
CTEQ66 Original asymmetry (MC input) 
Interation 1 
Interation 2 
Interation 3 



vtfon 



2 2.5 

absolute rapidity I y 



a) W asymmetry of MC input, pseudo-data and 
interation steps 
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LO MC input (CTEQ66 PDF, triangles) 
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compared a) to the reconstructed W 
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steps 1-5. b) depicts the ratio of the 
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iteration steps to the true W asymmetry 
of the pseudo-data, c) shows the ratio of 
the reconstructed W asymmetries of the 
iteration steps to the W asymmetry of the 
preceeding iteration step. 
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the first iteration). Except for the first iteration step, all of the other iteration steps agree 
very well with the reconstructed asymmetries from the preceeding iterations indicating the 
convergence of the iterative procedure. This demonstrates that the modified reweighting 
scheme is able to provide a good measurement of the W asymmetry. 

An indication for when the iteration should be stopped can also be taken from this 
last plot: As long as the ratio of the reconstructed W asymmetry of iteration step n to the 
reconstructed W asymmetry of iteration step n — 1 grows significantly closer to one with 
every iteration n, the iteration still converges. If the ratio does not change significantly 
from some constant numbers (as e.g. in the first bin of yw) or starts to deviate from one, 
the iteration starts to diverge and should be stopped well before this happens. 

All MC codes only simulate the physics process up to a certain order in QCD. The im- 
pact of higher order QCD effects was investigated by studying the analysis technique using 
NLO pseudo-data generated with Herwig++ and LO MC input generated with Pythia 
each using the same PDF set. Figure |l0| a) compares the true W asymmetries for the LO 
MSTW08 MC input and the NLO MSTW08 pseudo-data to the reconstructed W asym- 
metries. Again, the reconstructed W asymmetries agree well with the true W asymmetry 
of the pseudo-data, as also seen in figure [UJ b). The agreement is of the level of 5% for 
the first iteration. Statistical effects are larger due to the size of the data sets (1 million 
events for the NLO MSTW08 pseudo-data and 18 million events for the LO MSTW08 MC 
input). This can also be seen in figure |l(] c), where the deviations of the iteration step 
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from the previous iteration step are quite small for the first and the second iteration and 
then become systematically larger. Since this indicates a divergence, that grows with the 
number of iteration steps, the iteration should be stopped after the first or second itera- 
tion. Apart from the large statistical fluctuations, no significant systematic effects were 
encountered from higher order QCD. 

8. Comparison of Uncertainties on the Lepton and W asymmetry 

In order to assess the possible advantage of the direct reconstruction of the W asymmetry 
in comparison to a measurement of the lepton asymmetry, the expected statistical and 
systematic uncertainties on measurements of the two observables were compared to the 
uncertainty on the respective predictions arising from the PDFs. Here, for the W asymme- 
try, the new scheme as described in section ^ was used and acceptance cuts were applied. 
For the lepton asymmetry, the same acceptance cuts were applied as would be the case for 
the usual experimental determination of this quantity. 

Table [l] compares the relative uncertainties in detail for the W and lepton asymmetries 
including experimental uncertainties estimated using Gaussian resolution functions of 1% 
for the lepton PJp and 5% for the $t- 

Table [2] compares the uncertainties on the theoretical predictions from the PDFs (ctpdf) 
to the total uncertainties on the measurement (both with and without the experimental 
uncertainty due to the detector resolution). Here, the measured W asymmetry refers to 
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1.16 
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4.38 
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Table 1: Comparison of the relative uncertainties on the measured W asymmetry and the lepton 
asymmetry for a sample corresponding to an integrated luminosity of 268 pb -1 . The statistical 
uncertainty in the bins is denoted er s tat, uncertainties on acceptance corrections from PDFs are 
denoted <7 acc (c.f pi step 3). Uncertainties from detector resolution are denoted as tJdct- The total 
uncertainty neglecting <7d c t is labeled cr^ eas and the total uncertainty including resolution effects 



the W asymmetry extracted in the last step of the analysis technique (c.f section y, step 7) 
with the modified scheme described in section ||. Table ^ also gives the ratios of the errors 
on the measurement, <r^5r as , to opdf for the W asymmetry, 1Z(W), and lepton asymmetry, 
1Z(1). These ratios give an estimate of how well a measurement could be used to constrain 
the PDFs. Also shown is the ratio of these ratios, TZ{W/l) = 7Z(W)/7Z(l). This is an 
estimate of how the observables perform in comparison to each other, with a ratio of 0.3 
predicting the W asymmetry measurement to be only a third as well suited to constrain 
PDFs as the lepton asymmetry for all yw when the uncertainties without the experimental 
uncertainty are taken into account. The detailed numbers for different yw bins are given 
in table [l] and [2], in which they are calculated for 268.1 pb -1 integrated luminosity and are 
the relative uncertainties. 

The numbers are crude estimates since only the generator level is studied. However, 
the message from table |2| is clear: the measurement of the lepton asymmetry will be much 
more appropriate to constrain the PDFs, regardless of whether only the inherent theoretical 
uncertainties from the method are considered or whether experimental uncertainties are 
also taken into account. Measuring the W asymmetry has a power to constrain PDFs which 
is only 33% as good as a measurement of the lepton asymmetry, when only the inherent 
theoretical uncertainties on the method are taken into account (c.f. 7£(W / /7)<rpDF/c]^ eas )- 
Taking into account experimental resolution effects, the lepton asymmetry still performs 
twice as well compared to the W asymmetry measurement (c.f. lZ{W/l)a-pDY/o'^^ s )- The 
ratio IZw /Tt\ tends to decrease with higher luminosity. 

It should be noted again that the situation at the Tevatron is more favourable for the 
measurement of the W asymmetry. Since the lepton asymmetry at the Tevatron tends to 
be smaller than the W asymmetry, at the Tevatron the lepton asymmetry is favoured if 
the systematic uncertainties are negligible or of the same absolute size on the W and the 
lepton asymmetry. 



9. Conclusions 

In this paper the applicability of the analysis technique first developed for the Tevatron 
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Table 2: Comparison of the relative uncertainties on the measured W and lepton asymmetries and 
the uncertainties on the prediction arising from PDFs for a sample corresponding to an integrated 
luminosity of 268 pb -1 . The uncertainties on the theoretical predictions for the W and lepton 
asymmetry are denoted as cfpdf(W) and 0pdf(O respectively. Ratios of uncertainties for W and 
lepton asymmetries are denoted as R(W) or R(l) respectively. Finally double ratios 1Z(W)/1Z(1) 
are denoted as 1Z(W/1). 



by Bodek et al. [16|] to the LHC environment has been studied. The method is inherently 
less well suited for LHC conditions and a modification to the method has been proposed. 
This modification was found to be more robust and relies on MC input only for one step. 
The statistical and systematic uncertainties on the measurement of the W asymmetry 
were determined and compared to the PDF uncertainties on the prediction from theory. 
The expected experimental uncertainty on the measured W asymmetry is smaller than 
the theoretical uncertainty. However, the expected experimental uncertainty on the lepton 
asymmetry is even smaller in comparison to its theoretical uncertainties and thus better 
suited to constrain PDFs at the LHC. An estimate of the likely impact of experimental 
uncertainties due to the detector resolution were taken into account and the lepton asym- 
metry was still found to perform better than the W asymmetry by a factor of roughly 2. 
It is hence the opinion of the authors that early measurement at the LHC should focus on 
the lepton asymmetry rather than the direct W asymmetry to constrain the PDFs. The 
studies have been performed for a centre of mass energy of yfs = 14 TeV. However, the 
main arguments hold also at yfs = 7 TeV, so that these conclusions are not likely to change 
with the centre of mass energy. 
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